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TURBULENCE MODELING REQUIREMENTS, 
DEVELOPMENT PHILOSOPHY AND APPROACH 

• REQUIREMENTS 

• TURBULENCE MODELING IS A KEY ENABLING TECHNOLOGY 
FOR ALL PROPULSION RELATED CFD ACTIVITIES 

• FACTORS TO CONSIDER INCLUDE ACCURACY, CONSISTENCY, 
COMPUTATIONAL COST, AND EASE OF USE 

• TURBULENCE MODELS THAT CAN NOT BE INCLUDED IN 
PRODUCTION GRADE CFD CODES ARE OF LIMITED VALUE 
TO INDUSTRY 

• PHILOSOPHY 

• BASIC MODEL DEVELOPMENT IS BEST LEFT TO SPECIALIZED 
"CENTERS OF EXCELLENCE" 

• VARIOUS CLASSES OF MODELS NEED TO BE SUPPORTED SINCE 
NO SINGLE UNIVERSAL MODEL IS SHOWN TO EXIST 

• ESTABLISHING THE RANGE OF APPLICABILITY, ACCURACY, AND 
THE COMPUTATIONAL COST OF THE MODELS IS ESSENTIAL 


, TURBULENCE MODELING REQUIREMENTS, 
DEVELOPMENT PHILOSOPHY AND APPROACH (Cont.) 


• APPROACH 


• IDENTIFY KEY "CENTERS OF EXCELLENCE" AND ESTABLISH 
COLLABORATIVE RELATIONSHIP 

• ACQUIRE MODELS AND ASSESS PERFORMANCE FOR THE 
INTENDED CLASS OF APPLICATIONS 

• DELINEATE MODEL DEFICIENCIES AND INITIATE EFFORT TO 
REDUCE THEM 

• DEVELOP MODELS INTO STAND-ALONE MODULES 

• INCLUDE MODULES IN PRODUCTION CODES AND ESTABLISH 
BASELINE FOR APPLICATIONS 
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TWO MAJOR AREAS OF CONCENTRATION 


HIGH SPEED TURBULENCE MODELING (LEAD DR. DOUG LYNCH) 

• FOCUSED ON HIGH SPEED (M>1) PROPULSION (ROCKET 
AND AIRBREATHING) AND AERODYNAMICS 

• EMPHASIS ON 2-EQUATION PHENOMENOLOGICAL MODELS 
WITH NASA ARC AND LARC AS KEY TECHNOLOGY PARTNERS 

• LES WORK IN PLANNING STAGES WITH CTR 

LOW SPEED TURBULENCE MODELING (LEAD DR. ALI HADID) 

• FOCUSED ON LOW SPEED (M<1) AND ROTATING FLOW 
APPLICATIONS 

• EMPHASIS ON REYNOLDS STRESS PHENOMENOLOGICAL 
MODELS IN COLLABORATION WITH UMIST, ICOMP, CTR, AND UAH 

• LES WORK INITIATED WITH CTR 


HIGH SPEED TURBULENCE MODELING 

• EMPHASIS IS ON THE DEVELOPMENT OF ENGINEERING 
TURBULENCE MODELS FOR 

• HIGH SPEED AIRBREATHING PROPULSION SYSTEMS 

• THRUST CHAMBERS 

• VEHICLE AERODYNAMICS 

. APPROACH TAKEN IS BASED ON 2-EQUATION MODELS 

• DIFFERENT CLASSES OF 2-EQUATION MODELS STUDIED 

• k-e 

• k-<D 

• POINTWISER t 

• COMPRESSIBILITY EFFECTS AND TURBULENCE-CHEMISTRY 
INTERACTIONS MAJOR MODEL UPGRADE THRUSTS 

• COMPRESSIBILITY MODIFICATIONS FROM ARC 

• TURBULENCE-CHEMISTRY INTERACTION MODELS FROM LARC 

• USA AND GASP SERVE AS NUMERICAL PLATFORM 

• GASP - CHIEN, LAM-BREMHORST k-e, k-<o 

• USA - VARIETY OF k-e, k-a> 
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COMPRESSIBILITY EFFECTS 

• MIXING LAYER SPREADING REDUCED AT HIGH MACH NUMBERS 

• INCREASE DISSIPATION RATE OF k 

• DEFINE Ck 2 AS A FUNCTION OF TURBULENT MACH NUMBERa^HTT 

• ZEMAN MODIFICATION (1990) K 

• SARKAR (1990, 1991) AND WILCOX (1991) PROPOSALS 

• MODIFICATIONS OF ZEMAN AND SARKAR NOT RECOMMENDED 

• HEAT TRANSFER OVER PREDICTED NEAR SHOCK WAVES 

• LIMIT TURBULENT LENGTH SCALE L t TO MIN(k^ , J£y_\ 

(VUONG AND COAKLEY, 1987) £ C F3 /4 

• SEPARATION UNDERPREDICTED IN RAPID COMPRESSION OR 
STRAIN REGIONS 

• INCREASE acOR a SUNDER RAPID COMPRESSION (VUONG AND 
COAKLEY) 

• HEAT TRANSFER OVER PREDICTED FOR VERY COLD WALLS 
T v/T aw <0.1 (COAKLEY ) 

• CEBECI-SMITH ~ 60%, k-co ~ 40%, q-© ~ 10%, k-e ~ 30% 


TURBULENCE MODELS ADAPTED TO USA CODE 






TRANSITION MODEL 

COMPRESSIBILITY EFFECTS 


DAMPING 

BOUNDARY 

HIGH ORDER 


MIXING LAYER 

SEPARATION 

REATTACHMENT 

ALGEBRAIC 

m 

LOCAL 

CQHDflTONS 

POLYHQKBNAL 

AHHAL 

SPREADING 

EXTENT 

HEAT TRANSFER 

Baldwin- Lomax 

X 

3 Varslona 


X 

X 




fee 

1. Myong-Kawgl 

X 

X 

k*o 

X 


1. Sarkar (1991) 

2. Zaman (1990) 

4. YuongA 
Coaklay 

S. Vuong A 
Coaklay 






3. WUcox (1991) 

(1987) 

(1987) 

2. Chian (1982) 

X 

. X 

1:8 

X 


1-.2. 3. 

4. 

5. 

3. Jonas-Laundar 

X 

X 

kxO 

X 


1,2,3- 

4. 

5. 

(1972) 



d0 

k*0 
c * 0 
k»0 



1,2,3. 

4. 

S. 

4. Laundaf-Sharma 
(1974) 

X 

X 

X 


1,2,3. 

4. 

s. 

5. Huartg-Coaklay 
(1992) 

X 

X 

t» 2wjk,/y2 

X 







k.O 



1,2, 3. 

4. 

5. 

6. SpazJala-So-Zharvg 

X 

X 

i-2v,k,/y2 

X 





(1993) 








7. Lam-Bramhorst 

X 

X 

k«o 

CyaO 

X 


1,2,3. 

4. 

S. 

(1991) 









«. High fla 

X 

X 

Wall Function 

X 





La 









1. High fia Wilcox 
(I99ia) 

- 

- 

k*(J 
* i cki>| 

X 


1,2,3. 

4. 

5. 

2. Low R« Wilcox 

- 

- 

k*o „ 

»« 7.2\>j/y 2 

X 


1,2.. 3. 

4. 

5. 

(1991b) 








Q-<P 



kaO 

X 





Coaklay (1987) 



(DylO 






Ona-Equallon (Goldbarg, - 


kaO 






Two- Tima Seal# 1 992) 



<?) <K»tR,) 






On»-Equ*tlon Rj 


X 

= 0 X 






(Goldbarg 1993, 1994) dy dy 


109 



M = 9.2 PLAT PLATE FLOW 

. CHIEN k-e MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 


VELOCITY PROFILE 


TURBULENT VISCOSITY 



REF- G.T. COLEMAN AND J.L. STOLLERY, JFM 56: 741, "HEAT- TRANSFER FROM A 
HYPERSONIC TURBULENT FLOW AT A WEDGE COMPRESSION CORNER" 


MACH 7.05 PLOW OVER AXISYMMETRIC FLARE 

CHIEN k-(u MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 


AXISYMMETRIC FLARE 


WALL PRESSURE FOR AXISYMMETRIC FLARE 


5 = 2 cm 


M = 7.05 




D = 20 cm 


— — HO PREDICTIONS 
• ARC PREDICTIONS 
— — TEST DATA 



REF- M I KUSSOY AND C.C. HORSTMAN, "DOCUMENTATION OF TWO- AND 
' THREE-DIMENSIONAL HYPERSONIC SHOCK-WAVE TURBULENT 
BOUNDARY LAYER INTERACTION FLOW," NASATM 1-01075. 







MACH 7.05 FLOW OVER AXISYMMETRIC FLARE 

CHIEN k-co MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
COMPRESSIBILITY MODIFICATIONS 


WALL HEAT TRANSFER FOR AXISYMMETRIC FLARE 



x (in) 


MACH 8.6 FLOW OVER COLD WALL WEDGE 


M = 8.6 

Tw/Taw = 0.065 

THREE STUDIES 

1 . CHIEN k-e MODEL WITH RAPID COMPRESSION AND LENGTH SCALE 
CORRECTIONS AND WITH AND WITHOUT MIXING LAYER TREATMENT 

2. HIGH-Re k-co MODEL WITH VARIOUS AIR CHEMISTRY MODELS 

3. BALDWIN-LOMAX TURBULENCE MODEL USING WALL AND LOCAL 
DAMPING 
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MACH 8.6 FLOW OVER COLD WALL WEDGE 

CHIEN k-e MODEL WITH AND WITHOUT MIXING LAYER TREATMENT 
HEAT FLUX CALCULATIONS 



MACH 8.6 FLOW OVER COLD WALL WEDGE 

HIGH-Re k-0) MODEL WITH VARIOUS AIR CHEMISTRY MODELS 


HEAT FLUX CALCULATIONS 
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MACH 8.6 FLOW OVER COLD WALL WEDGE 

BALDWIN LOMAX, k-e, k-co MODEL COMPARISONS 


HEAT FLUX CALCULATIONS 



LOW SPEED TURBULENCE MODELING 


• EMPHASIS IS ON THE DEVELOPMENT OF ENGINEERING 

TURBULENCE MODELS FOR 

• ROTATING MACHINERY 

• FLOW IN DUCTS AND MANIFOLDS 

• REACTING FLOWS 

• APPROACH TAKEN IS TO 

1 . SYSTEMATICALLY ASSESS EXISTING PHENOMENOLOGICAL 
MODELS USING COMMON NAVIER-STOKES SOLVER 

2. IDENTIFY, DEVELOP AND VALIDATE MODEL UPGRADES 
COMMENSURATE WITH OBSERVED FLOWPHYSICS 

3. DEVELOP SELF-CONTAINED TURBULENCE MODEL DECKS 
(MODULES) THAT CAN BE INTEGRATED WITH NAVIER-STOKES 
SOLVERS 

4. PROVIDE GUIDANCE TO EXPERIMENTAL AND THEORETICAL 
RESEARCH IN TURBULENCE MODELING FOR ENGINEERING 
APPLICATIONS 
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TURBULENCE MODELS BEING ASSESSED 


PHENOMENOLOGICAL SINGLE POINT CLOSURE MODELS 


2-EGUATION MODELS 
K- e (SKEM) 


SINGLE-SCALE 


MULTI-SCALE 


2-EQUATION MODELS 
K- e (MKEM) 


ALGEBRAIC REYNOLDS 

STRESS MODELS STRESS MODELS 
(ASM) (RSM) 


NEAR-WALL TREATMENTS INCLUDE (WHERE APPROPRIATE) WALL 
FUNCTIONS, MULTILAYER MODELS, AND LOW-REYNOLDS NUMBER 
APPROXIMATIONS 


TURBULENCE MODEL DECK STRUCTURE AND 
INTEGRATION WITH NAVIER-STOKES SOLVER 
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PROJECT WELL UNDERWAY 


• TEAM 

• MODELS PROVIDED BY UMIST, LERC/ICOMP, ARC/CTR 

• MODULE DEVELOPMENT BY ROCKETDYNE 

• MODULE TESTING BY ROCKETDYNE (REACT, USA) AND 
UAH (MAST) 

• MODEL UPGRADES BY ROCKETDYNE, UMIST, ARC/CTR 

• APPLICATION BY ROCKETDYNE TO TURBOPUMP 
COMPONENT (E.G. IMPELLER) ANALYSIS 

• 2-D MODULES COMPLETED, TESTED, AND RELEASED 

• SINGLE SCALE k-e 

• MULTI SCALE k-e 

• ASM 

• RSM 


• 3-D MODULE DEVELOPMENT IN PROGRESS 


NONLINEAR ALGEBRAIC-STRESS MODEL 

VORTEX SHEDDING FROM RECTANGULAR CYLINDERS (DURAO, et al) 

PARTICLE STREAKLINES 



MEAN AXIAL VELOCITY ALONG CENTERLINE MEAN KINETIC ENERGY ALONG CENTERLINE 
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ROTATION MODIFIED k-e MODEL 

BACKWARD FACING STEP (DRIVER AND SEEGMILLER) 


STREAMLINE CONTOURS 


0.00 



MEAN AXIAL VELOCITY AT X/M=4 


RADIAL TURBULENT INTENSITY (w/i^ref) 




ALGEBRAIC STRESS MODEL 

CONFINED COAXIAL SWIRLING JET FLOW (ROBACK AND JOHNSON) 


GEOMETRY 


STREAMLINE CONTOURS 




REYNOLDS STRESS MODEL (LRR - MODEL) 

BACKWARD FACING STEP (DRIVER AND SEEGMILLER) 

STREAMLINE CONTOURS 



MEAN AXIAL VELOCITY AT X/H=4 AXIAL TURBULENT INTENSITY AT X/H=4 



CONCLUDING REMARKS 

• PROGRAMS (BOTH COMMERCIAL AND GOVERNMENT) EMPLOY 
NEW TECHNOLOGY ONLY WHEN IT PROVIDES "ADDED VALUE" 

• REDUCED DEVELOPMENT COST 

• INCREASED RELIABILITY AND PERFORMANCE 

• ENHANCED MANUFACTURABILITY 

• THE NEW TECHNOLOGY WE OFFER IS THE COMPUTATIONAL 
ENGINEERING TOOLS FOR PRODUCT DESIGN AND ANALYSIS 

• THESE TOOLS ARE THE END PRODUCT FOR ALL ENABLING 
TECHNOLOGY DEVELOPMENT 

• PRE- AND POST PROCESSING 

• ALGORITHMS AND NUMERICAL PLATFORMS 

• PHYSICAL MODELS (E.G. TURBULENCE AND CHEMISTRY) 

• FAILURE OF ANY ENABLING TECHNOLOGY JEOPARDIZES THE 
PERFORMANCE (VALUE) OF THE TOOL 

NOW MORE THAN EVER, THERE IS A NEED FOR CLOSER 
COLLABORATION AND COOPERATION BETWEEN GOVERNMENT, 
INDUSTRY, AND RESEARCH INSTITUTIONS TO ENSURE 
MAINTENANCE OF COUNTRY'S TECHNOLOGY BASE 
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